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We study a one-dimensional system that consists of an electron gas coupled to a spin-1/2 chain 
by Kondo interaction away from half-filling. We show that zero-temperature transitions between 
phases with "small" and "large" Fermi momenta can be continuous. Such a continuous but Fermi- 
momentum-changing transition arises in the presence of spin anisotropy, from a Luttinger liquid 
with a small Fermi momentum to a Kondo-dimer phase with a large Fermi momentum. We have 
also added a frustrating next-nearest-neighbor interaction in the spin chain to show the possibility of 
a similar Fermi-momentum-changing transition, between the Kondo phase and a spin-Peierls phase, 
in the spin isotropic case. This transition, however, appears to involve a region in which the two 
phases coexist. 

PACS numbers: 75.30.Mb, 71.10.Hf, 71.10.Pm, 71.30.-|-h 



I. INTRODUCTION 

One of the key issues in the theory of heavy fermion 
materials is the size of the Fermi surface. 

It is traditional to think that each spin-1/2 localized 
moment, through the Kondo effect, contributes a Kondo 
resonance. These Kondo resonances, which are fermionic, 
combine with the conduction electrons to form (very 
heavy) quasiparticles. In this way, the local moments 
ultimately become a part of the low-energy electronic 
fluidi This is reflected in the size of the Fermi surface: 
the Fermi momentum kp is "large" in the sense that the 
Fermi surface encloses a volume that counts the number 
of both conduction electrons and local momentsi^ This 
is in contrast to a "small" Fermi momentum kp, which 
would correspond to a Fermi surface that has a volume 
counting only the number of conduction electrons. While 
some questions have been occasionally raised, including 
the exhaustion problem noted by NoziereS)^ this pic- 
ture has been thought to apply everywhere in the heavy 
fermion phase diagram i^i^ 

Recently, the question about the size of the Fermi sur- 
face has gained renewed importance in the light of the 
anomalous behavior observed in the vicinity of the an- 
tiferromagnetic quantum critical point (QCP)ji Particu- 
larly noteworthy are the experiments in CeCug-xAux, 
YbRh2Si2, and related systemSj*'^'^° which observed 
non-Gaussian behavior. These experiments have in- 
spired the development of locally critical quantum phase 
transitions^^ and related theoretical pictures.^^ In these 
pictures, the paramagnetic heavy fermion metal phase, as 
usual, contains fully developed Kondo resonances, with a 
large Fermi momentum kp. It is argued that the Kondo 
effect can be destroyed as the system approaches a QCP 
and goes into an antiferromagnetic metal phase, by some 
competing processes: the local moments in such a lat- 
tice setting are also coupled to some fluctuating magnetic 
field, which impedes the Kondo effect. Once the Kondo 
effect is destroyed, the local moments are no longer a 
part of the electron fluid and the Fermi surface becomes 



small. Other theoretical works have proposed param- 
agnetic phases in which the Kondo effect is destroyed 
due to the fractionalization taking place in the localized- 
spin component. One of the important questions that 
these theoretical developments raised is whether quan- 
tum transitions between phases with very different Fermi 
momenta is necessarily first order or can be continu- 
ous. Only in the latter case can transitions between such 
phases serve as a mechanism for the non-Fermi liquid 
behavior observed in the heavy fermion metals. 

In this paper we address this issue in models in one di- 
mension, where they are amenable to controlled compu- 
tation. We show that continuous transitions from small 
to large Fermi momenta can indeed take place in a one 
dimensional Kondo lattice model away from half-filling. 

We now discuss a few additional background issues be- 
fore we go into the details of our analysis. 



A. Small and large Fermi momenta 

We prefer to speak about small and large Fermi mo- 
menta as opposed to small and large Fermi volumes. 

The situation in the paramagnetic phase is straightfor- 
ward and is standard in the heavy fermion literaturei^ 
Whenever the Kondo effect takes place, the localized 
spins are part of the low-energy electron fluid. The Fermi 
momentum in this case is the same as that when each lo- 
calized spin is replaced by an electron orbital and, in ad- 
dition, with all two-particle interactions turned off (i.e., 
a non-interacting Anderson lattice model). We call this 
Fermi momentum a large Fermi momentum. On the 
other hand, in the absence of the Kondo effect, local- 
ized spins — being charge neutral — are not part of the 
electron fluid. The Fermi momentum is that of the con- 
duction electrons alone with all interactions turned off. 
This is defined as a small Fermi momentum. Clearly the 
Fermi volume (more precisely, Fermi length in one di- 
mension and Fermi area in two dimensions) is different 
in the two cases. A transition between a large-Fermi- 
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momentum phase and a small-Fermi-momentum phase 
involves a jump in Fermi volumes. (This is distinct from 
other electronic topological transitionsii^ For example, 
in the case of Lifshitz transitions in metals, the Fermi 
volume changes continuously, while the behavior of the 
density of states is anomalous.) 

The situation is more subtle in an ordered phase that 
breaks translational invariance. This includes states with 
either an antiferromagnetic order or a spin-Peierls order. 
In the case of even-integer-fold breaking of discrete trans- 
lational invariance of the Hamiltonian, the Fermi volumes 
(Fermi area in two dimensions and Fermi length in one 
dimension) in the reduced Brillouin zone would be the 
same regardless of whether the local moments are a part 
of the electron fluid or not. 

However, for commensurate ordering, large Fermi mo- 
mentum and small Fermi momentum can still be defined 
as follows. The key distinction is again whether local 
moments are a part of the electron fluid. We define a 
large Fermi momentum to be one that is associated with 
each localized spin being replaced by an electron orbital, 
all two-particle interactions turned off, but in the pres- 
ence of an infinitesimal static field that characterizes the 
symmetry breaking. (In a commensurate antiferromag- 
netic metal phase, this field would be a static staggered 
magnetic field.) Likewise, we define a small Fermi mo- 
mentum to be that associated with conduction electrons 
alone, with all two-particle interactions turned off, and in 
the presence of the same infinitesimal static field. In di- 
mensions higher than one, the corresponding large Fermi 
surface and small Fermi surface have different topolo- 
gies and, hence, represent distinct quantum phases. In 
one dimension, on the other hand, the only distinction 
between the two is the actual location of the Fermi mo- 
mentum; the two phases may not necessarily be distinct 
if there are no additional order parameters that differen- 
tiate them. 



B. The one dimensional Kondo lattice 

The system of interest consists of a one-dimensional 
electron gas (IDEG) coupled to a one-dimensional spin- 
1/2 chain via Kondo interaction. The spins in the 
spin chain Tj interact with each other by nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) inter- 
actions. The IDEG is characterized by Fermi energy ep 
and Fermi momentum kp. We take the lattice constants 
for both the conduction electrons and spin chains to be 
equal to a. Away from half- filling, kp ^ 7r/2a so that the 
IDEG and the spin chain are incommensurate with each 
other. 

The nature of Fermi momentum in various phases of 
such a one-dimensional Kondo lattice has been addressed 
in the past. That the Fermi momentum in the Kondo 
phase is indeed kp is supported by a generalized Lut- 
tinger's theorem, derived in a way analogous to that 
for the Lieb-Schultz-Mattis theorem, and by the numeri- 



cal calculations--^'- However, indications for phases with 
small Fermi momentum have been shown in more recent 
density matrix renormalization group calculations^ and 
the effects of spin dimerization in this context have been 
recently discussed. 

In absence of Kondo interaction, the IDEG and the 
spin chain are completely decoupled. In order to study 
the low-energy excitations, one can take the continuum 
limit so that IDEG will be described as a Luttinger liq- 
uid. 

For large NN interaction, the spin chain will also be- 
come a chargeless Luttinger liquid, with NNN interac- 
tion serving as a perturbation. When the perturbation 
becomes sufficiently large, the system undergoes a phase 
transition into the spin-Peierls phase. When the NNN 
interaction becomes much larger than the NN interaction, 
the spin-Peierls ordering persists, but it becomes more 
convenient to describe the system in terms of a "zigzag" 
modelf22iS in which the odd spins and the even spins 
form two separate chargeless Luttinger liquids and NN 
interaction introduces a weak coupling between them. 

The spin-Peierls phase is characterized by staggered 
bond order parameter, (r, • Tj+i) ^ (~1)"'- Thus, the 
spins form dimers, the translational symmetry is broken, 
and in certain limits one can even represent the ground 
state as a chain of singlets formed by nearest-neighbor 
spins. This phase has an exponentially small spin gap. 
The IDEG remains a Luttinger liquid with Fermi mo- 
mentum kp. 

Non-vanishing Kondo interaction introduces a number 
of novel phases in the systemfi2iSiS£iS& however, we will 
be primarily interested in limit where the Kondo inter- 
action is weak. In this case and for an incommensurate 
lattice, a relevant Kondo coupling leads to a phase^^ in 
which the electrons and some of the spins form Kondo 
singlets s{xj) ■ tj, where s (x) — (x) crtp {x) and 
tp (x) , tp^ {x) are the annihilation and the creation op- 
erators in the IDEG. The Kondo singlets are further 
dimerized so that the phase is characterized by an or- 
der parameter <I>k. 

Due to the coupling to the charge sector, this order is 
associated with gapless charge-density modes O at 2kp, 

0{X^) ^ S{X^) ■ Tj, (1) 

so that the expectation value of this operator varies as 
(O) ~ g^'^kpXj ^ where kp ^ kp + 7r/2a [a more detailed 
expression will be given in Eqs. (|20I21() later in the text]. 
For these modes the charge and spin quantum numbers 
are both zero. Note that the charge-density modes are 
described by a local operator, however, due to the fact 
that they embed the Kondo singlets, the order parame- 
ter <I>K becomes nonlocal. As a result, the translational 
symmetry of the system is not broken and the order is 
hidden. Furthermore, there is an exponentially small spin 
gapf2i2& which exists due to the fact that the order pa- 
rameter is nonlocal and one cannot distinguish paired 
spins from unpaired ones. The existence of the gap- 
less charge-density modes at 2kp implies that the system 
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should be described as one with a large Fermi momen- 
tum. 

We will argue in the following for a new co-existence 
phase. In the Kondo dimer phase the spins can further 
dimerize, producing a mixed phase in which the Kondo 
dimer order coexists with spin-Peierls order. This phase 
has both a broken translational symmetry and a large 
Fermi momentum. 

In addition, we will show that, away from the spin- 
isotropic case, the decoupled Luttinger liquid phase has 
a finite range of stability. This phase is paramagnetic 
and has a small Fermi momentum kp. 

We now proceed to the quantitative analysis of the 
model, which appears to possess all the mentioned 
phases. 



II. MODEL 



The Hamiltonian of the system is 



We will study this model in the limit Jk, J2 ^ Ji <C 
e F . Then we can take a continuum limit for the spin chain 
and describe it also as a one-dimensional free chargeless 
fermion system with Fermi energy of order of Ji. Hence, 
we can perform a similar decomposition of Tj into right- 
and left-directed currents and 3'^, as well as the stag- 
gered spin ny = ( — 1)-' n"^: 



(7) 



W — WlDEG + Wii 



Tin 



(2) 



Likewise, in the limit Ji <C J2, the continuum limit of 
the spin chain can be taken using a system of two mu- 
tually noninteracting fermion gases characterized by the 
Fermi energies of order of J2 (the "zigzag" model), with 
Ji serving as a perturbation. ^^-^^ 

For the single impurity Kondo problem and for the 
commensurate Kondo lattice problem, the most relevant 
part of Kondo interaction is nj • n*. However, for the in- 
commensurate case — which is what happens away from 
half-filling — this interaction becomes irrelevant due to 
the oscillating factors, and the only component that is 
not irrelevant is the forward-scattering one. 



where TYideg is the kinetic term of non-interacting one- 
dimensional electron gas (IDEG) with Fermi energy ep, 
Tinn is Heisenberg nearest-neighbor interaction between 
the spins in the one-dimensional 1 /2-spin chain. 



(3) 



Tinnn is Hciscnberg next-nearest-neighbor interaction in 
the spin chain. 



Ti-n 



T7+2, 



(4) 



'^anis introduces anisotropy in the spin chain, whose ex- 
act form is not essential provided that the spin chain by 
itself {Hnn + 'Wnnn + '^anis) is gaplcss. We cau, for in- 
stance, take Hanis = "-^'^I'^/'^Z+n ^i^i^ < A < 1, which 
is an easy-plane (XY) anisotropy. [For a discussion of 
the role of anisotropy, see Appendix ^] Finally, Hk is 
Kondo interaction between the IDEG and the spin chain, 



(5) 



f - jf \^ 



T-/^ — J 



[3% +31) {J], + 31), 



(8) 



which is marginally relevant in the SU (2)-symmetric 
case. 



Our next step is Abelian bosonization 
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which we represent the fermionic fields of IDEG as 
Ra (x) — [Fa/ V2Tra) exp [9a — (j)cr]} and Lcr (x) = 

(Fa/V^TTa) exp{-iV7r [6*0- + (/)„]}, a =t,i, where (j)^ (x) 
are the bosonic fields, 9a (x) = Ha {x') dx' are the 
fields dual to 4>a, H^ (x) are the conjugate fields satis- 
fying [(jja [x) ,Iia' [x')] — —i5aa'5{x — x'), and Fa are 
the Klein factors satisfyingSL^^ F^Fa = FaF^ = 1 and, 
for a + a', F]F,, = -Fa>Fl and = -Fa'Fa. 

The bosonic fields are further re-expressed in terms of 
the spin fields [x) = ((/)| — <^{) /\/2 and charge fields 
(j)'^ [x) = ((/)| -I- (j)i) / v^, with similar expressions for their 
duals 0* and 9'^. The spin chain is bosonized in a similar 
manner and we will denote its spin fields as (jf (x). 

The charge sector of the model separates from the 
spin one. Since localized spins do not carry charge, only 
IDEG contributes to this sector, which is described by a 
Gaussian model. 



We will focus on the antiferromagnetic case {Jk > 0). As 
a reminder, here Tj are the spins in the spin chain and 
s (x) = {x) fip (x) is the spin density of the electron 
gas at X. We decompose tp (x) into right- and left-moving 
components, ipa (x) = Ra {x) e^'^"^ + La {x) e"**-'^^, a = 
t , i , so that the spin density becomes 



s{x) 



(6) 



where J 



R 



and n'' 



1 



dx 



{dr^^)' + vl{d^4>^y 



(9) 



Therefore, we will focus on the spin sector. We will 
study the system in the general case when the spin chain 
can be anisotropic. Although we assume that the mi- 
croscopic Kondo interaction is isotropic (we will com- 
ment on the more general picture later), it is necessary 
to keep track of its diagonal component separately 
from the orthogonal component jj^ . As wc will show, in 
the process of the renormalization-group (RG) flow, the 
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originally isotropic Kondo interaction usually becomes 
anisotropic and restores the isotropy only near the fixed 
point of the Kondo phase. 

The resulting Lagrangian becomes 



K ' 



(10) 



where 
Co 

Csp — 



dx 



1 



2KrVr 

Jsp 



(tto) 



dx 



dx cos 



Ztt 



2J 



y dxcos \^p2^<if^ cos (^\/2^<?!)^ 
2^(r - 



(Tra) 
X cos 



(11a) 
(lib) 

(11c) 
(lid) 



Vs and Vr are the Fermi velocities of the IDEG and the 
spin chain, respectively, a is the lattice constant, and 
Jsp is the backscattering interaction, which favors spin- 
Peierls phase when it is positivejSi For the spin-isotropic 
case (when A = 0), Vr = 7rJi/2 Jsp = ! {J-i — J^) — 
Ci(J2- J|), where Ci ~ 1.72 tt^ and ~ 0.24 Ji,S 
Again, this assumes that J2 ^ Ji . If we increase J2 to 
the point when Jk, Ji <^ J2, it will be more appropriate 
to describe the system using the "zigzag" representation 
so that for both spin sub-chains v-r ^ J2 and Jsp ''^ Ji- 
In addition, finite A modifies the expressions for both v-r 
and J sp- 
in general, the value of Kr reflects the fact that the 
spin chain is anisotropic [we will assume that it has XY 
anisotropy and will generalize to Ising anisotropy in Ap- 
pendix 0. This substantially affects the scaling dimen- 
sions of the terms in Eq. H10|) . in particular, the backscat- 
tering and transverse Kondo interactions are no longer 
marginal. 

The Kondo part of this Lagrangian contains a term 
cos \/27r ((/)* — (j)^) cos \/27r (0* — 6''^), which is strongly ir- 
relevant. Thus, we can omit it in , leading to: 



J: 



/-L 
K 



{■nay 



dx cos V2tt (6^ -f ( 



I cos 



(12) 

This expression is quite remarkable. We observe that 
exchanging the bosonic fields with their duals ^ 
0*, 0"^ —9'^ and simultaneously replacing Ks^t with 
leaves the Lagrangian £0 + invariant, but 
induces changes in and £sp- This should manifest 
itself not only in the RG equations, but also in the low- 
energy physics. 

Although interaction jj^ alone is marginal and J^^ 
alone is marginal or irrelevant, these constants effectively 



modify the scaling dimensions of each other so that they 
both can become relevant. 



III. RENORMALIZATION-GROUP EQUATIONS 

In order to treat the combined effects of the various 
interaction terms appearing in the model Eqs. pOlllfl . 
we carry out an RG analysis. We will use a Coulomb-gas 
expansion)2^i2^ in which the interactions associated with 
the cosine terms {J^j^ and Jsp) appear as fugacities, and 
those associated with the quadratic terms {J^, Kr and 
Kg) specify stiffness constants. The quadratic parts of 
the Lagrangian Eqs. (|llalllcp can be diagonalized. 

The general aspects of the derivation of the RG equa- 
tions are standard^^i^ and will not be detailed here. 
What is nontrivial here, however, is the fact that we 
have two kinds of fugacities (those associated with J^"^ 
and Jsp), and the way they are coupled with each other 
is the key to the determination of the phase diagram. 
Below, we document some of the details regarding this 
coupling. 

The RG equations are found to be 



dl \ 2 2 2Ks 2Kr ' 



uKsKryi_Vz 



^Vz 2 



dysp 
dl 

dKr 
"dT 

dKs 
dl 



^2{l^Kr)ysp, 

= (1 - yl - 2Klylp, 



(1 - Kl) yl 



(13a) 
(13b) 

(13c) 

(13d) 

(13e) 



and we have introduced the 



where u = \/2ttVs/ (vg 4 

dimensionless coupling constants y± = Jk' / (St) 



y. - 47 (2^) 



3/2 , 



3/2, 



^3/2, 



andysp = Jsp/47r'^"wT-. Note that 
the coefficients of the y^ and y^p terms in Eqs. H13dll3e() 
actually depend on the method of regularization, how- 
ever, those of the y^ term must preserve the invariance 
of the equations (for ysp = 2/2=0) with respect to the 
change of variables Ks — > l/Kg and Kr — > 1/Kr due 
to the symmetry of Eq. 112|l discussed above. At the 
spin-chain isotropic point with ysp = 0, the flow is along 
I'Cg — Kr — 1. Far away from this point the Kondo in- 
teraction is irrelevant and its flow is dominated by the 
first term in Eq. (|13a|) . 

There is no relevant correction due to yz or y± on the 
right-hand side of Eq. (|13c|) . In the third order of per- 
turbation we have found the following term there: 



-27r 



vi {v, + 2vr) K,K^ 2 



Vr {Vs + Vr) 



-VzVSP- 



(14) 
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A similar expression has been obtained for a different 
spin-Peierls order as welli^ The fact that this correc- 
tion makes the backscattering interaction more irrelevant 
means that there is a competition between the Kondo- 
dimer and spin-Peierls phases and not a mutual enhance- 
ment. If there was enhancement (as in the case of melting 
in two dimensions), the two primary phases would be one 
with both kinds of order being present and one with no 
order. In our case, the mixed phase can only be inter- 
mediate, while the pure Kondo and spin-Peierls phases 
play major roles. Remarkably, in the second order of 
perturbation even the direct competition vanishes. 

In absence of either the backscattering Jsp or 
transverse Kondo J^"*" interactions, the remaining z- 
component of Kondo interaction is strictly marginal. 
However, the theory is still Gaussian21 and, by using a 
linear transformation of fields and (j)'^, we have ob- 
tained (see Appendix 0) an exact solution describing 
two decoupled liquids characterized by new spin wave 

velocities v'^, v'^, with corrections of order of (^Jji^ for 

J 1271 <^ Vs,Vt- Thus, we still have a Luttingcr-liquid 
state with four "Fermi points" in the spin sector, two 
at ztkp and two at ±7r/2a, which means that small J^? 
does not change the behavior of the system qualitatively. 
In each of the decoupled liquids the particles are par- 
tially in the spin sector of the original electron liquid and 
partially in the original localized-spin liquid. Increasing 
enhances this mixing until one encounters a singu- 
larity, which could be a signature of a transition into 
the "Toulouse point phase" with Fermi momentum kp 
or the development of quasi-long-range spin-density-wave 
order. '^^ 

There is a process that is not being taken into account 
in Eqs. ifT^ . During the RG flow, a term dxO^dxd'^ is 
generated. This term has a physical meaning of an in- 
teraction between the z-components of physical spin cur- 
rents of the electrons and the localized moments. Its 
coupling constant (let us denote it y^) renormalizes as in 
Eq. (|13b|) . with the replacement of with K:^^, and it 
makes a contribution to f3{yi_) similar to one that has 
already been made by y^. As long as y_L, ?/z ^ 1, this 
correction only negligibly affects the flow of y±. However, 
in the Kondo phase all three coupling constants become 
approximately equal to each other and renormalize to 
infinity. Taking y^ into account rigorously might affect 
certain numerical factors in the critical exponents, but 
it should not change the critical behavior qualitatively, 
therefore it has been neglected in Eqs. Ijl^l) . 

Let us now focus on the equations Eqs. H13all3bll3d|) 
involving only Kondo interaction in presence of 
anisotropic spin chain (which also corresponds to the case 
when ysp — and Ks = 1). We aim at determining 
the separatrix which divides the region where the Kondo 
coupling {y±) is irrelevant from the region where it is rel- 
evant. For small values of the bare Kondo coupling, we 
will see that the anisotropy x = Kr — 1 is small on the 
separatrix. For small a;, the RG equations simplify to the 



following: 



dy± 
dl 
dyz 
dl 
dx 

H 



1 



2 

uy'i 



X y±+ uy±yz, 



-2xyl. 



These equations have the following first integrals: 



(15a) 
(15b) 
(15c) 

(16a) 
(16b) 



The physical meaning of the constant Ci can be derived 
from the initial condition that the Kondo interaction is 
isotropic at the beginning of the flow, when x = Xq and 
y± = yz ~ yo, thus, Ci = Xq/4. While C2 can be un- 
derstood as the value oi x = x such that yz{x) = 0, 
it is only a formal definition, since yz never vanishes or 
changes sign for the mentioned initial conditions. There- 
fore, C2 = X actually determines the value of yo. 

The equations Eq. I|16|) can be written in a more con- 
venient form, 



± V: 



X 

T 



1/2 



yz 



U X 

— in — , 

2 X 



(17a) 
(17b) 



which represents the lines of the flow of the RG charges 
2/_L and yz (Fig.P). As we see, dy±/dl (and consequently, 
dy±/dx) vanishes at x = x, where 



1/2 



(18) 



For xo <C u, one can approximate x :^ xo — Xq/2u^. 
There is also a flow line that separates the lines ending 
at y_L = from the lines flowing towards y± = 00. This 
separatrix is determined by 



1 

a;exp 1 — 
u 



(19) 



which is aiscp — xq + Xq/2v? for xq <^ u. 

It follows that in general, even if initially y±^ = yz ~ yo, 
the Kondo interaction becomes anisotropic as a result of 
the RG flow such that y± < yz- In particular, at the 
line of Luttinger- liquid fixed points x < x < xq, y± — 0, 
the diagonal component of Kondo interaction yz is finite. 
However, if the lines flow towards the fixed point of the 
Kondo phase x ^ 0, y± =00, the isotropy is eventually 
restored, yz ^ y±- 

A nontrivial consequence of the presented calculation 
is that the flow along the separatrix ends at a point 
where not only Kondo interaction is anisotropic, but also 
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FIG. 1: The RG flow of J^^ vs. x = Kr - I for Jsp = 0. 
At xo the Kondo interaction is isotropic, which is the initial 
condition. 
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FIG. 2: The RG flow of Jsp for initially smaU Jk- 
simultaneous flow along Jk is not shown. 



The 



spin-chain XY anisotropy is characterized by a nonzero 
value X. This value is nonuniversal and is related to the 
strength of the z-component of Kondo interaction at the 
point. These features are unique to our problem and have 
to do with the fact that both and Kr affect the condi- 
tions for the development of the Kondo phase. However, 
since this transition involves the change of the behavior 
of one fugacity, we expect that the universal properties 
will be of the Kosterlitz-Thouless type. Indeed, by in- 
troducing e — ^x^ — uyz, the RG equations IjlSI) can be 
expressed in the standard Kosterlitz-Thouless form. 



IV. PHASES AND FERMI MOMENTA 

First, we consider the quantum phases as characterized 
by the fixed points of the RG flow. 



irrelevant in the RG sense. While the details of this mech- 
anism differ from those of either Ref. 1 1 (fluctuating mag- 
netic field) or Rcfs. (gapped spin liquid), the effect is 
similar: they lead to the destruction of the Kondo effect. 

The only subtlety in the calculation of correlation func- 
tions in this phase is the J^^-coupling, which is marginal. 
It mixes the spin degrees of freedom from the spin chain 
and from IDEG. This mixing can be straightforwardly 
treated by introducing a new basis, leading to two de- 
coupled spin branches with renormalized velocities. The 
details are given in Appendix^ The spin excitations are 
gapless at both 2kF and 7r/a, which is due to the fact that 

f Z 

does not induce a magnetic order and, consequently, 
does not create new spin excitations. 



A. Decoupled Luttinger liquid with small Fermi 
momentum 

For spin-isotropic systems, the decoupled Luttinger liq- 
uid phase is unstable. The situation is drastically differ- 
ent when we go to the anisotropic case. For simplicity, we 
consider Kg = 1 corresponding to the case in which the 
ID electron gas is non-interacting. If initially ^ ysp 
and Kr 3> 1, the backscattering interaction ysp rapidly 
renormalizes to zero (Fig.EJ. After that we only need to 
consider the flow of the yx- For Kr > 1, there is a finite 
range of the Kondo interactions (both and over 
which ypc renormalizes to zero. 

The impedance to the Kondo effect here comes from 
the formation of a Luttinger liquid among the localized 
spins. The latter changes the scaling dimension of the 
Kondo interaction and, over an appropriate region of the 
interaction parameter space, renders the Kondo coupling 



The charge sector is completely independent of the 
J^^-coupling, therefore, the number of particles in this 
sector does not change at all. Consequently, the Fermi 
momentum is not affected and the charge excitations are 
gapless only at 2kp. In addition, the single-electron ex- 
citations are gapless only a,t kp: while changes the 
shape of the spectral function, it does not induce gapless 
single-electron excitations at other wavevectors. In this 
sense, the Fermi momentum is small and is equal to kp. 

We should stress that there is no violation of the "gen- 
eralized Luttinger theorem" as specified in Ref. 13. This 
theorem states that the twist operation (U) will introduce 
a gapless state, U\gs), whose momentum differs from that 
of the ground state |gs) by 2kF + tt/q. In the decoupled 
Luttinger liquid, a part of the momentum (2kp) mea- 
sures the change to the IDEG and the other part (tt/o) 
measures the change to the spin chain. 
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B. The Kondo phase with large Fermi momentum 

Another RG fixed point is given by = cxd. As dis- 
cussed earfier,^^ this corresponds to the spin-gap phase 
with ^ {(j)" + (j)'') = (2ni-hl)7r, ^/2^ {6" - 6'') 
27rn2, where ni_2 are integers, or with A/27r (0" 4- (jf) — 
27mi, ^/2^T{e^ - 6^) = (2^2 + 1)tt. This phase is pre- 
cisely the state in which Kondo singlets n'' • nj form an 
order parameter 



$1. 



cos V 27r 



cos V 27r 



(20) 



This order parameter is associated with the gapless 
charge-density mode Eq. at 2k*p^ 



O {x) - $Ke' 



(21) 



In the left-hand side of Eq. , one can replace O {xj ) 
with —O (xj) O (xj+i), which follows from J^~^n^~ ^ n^^ 
and similar relations. Physically this new formula means 
that the spin sector of the Kondo phase contains the 
dimers of the Kondo singlets. However, the order pa- 
rameter $K also includes a contribution from the charge 
sector exp (— i-\/27r0^) . 

Only relative values of the phase bosonic fields at- 
tain nonvanishing expectation values, while the fields 
themselves remain fluctuating. This picture does not 
depend on whether the original Kondo interaction was 
anisotropic or not, thus, Kondo dimer phase restores the 
isotropy of the Kondo interaction. 



C. The spin-Peierls phase with a small Fermi 
momentum 

Yet another RG fixed point is given hy ysp — +oo and 
yK — 0, corresponding to the spin-Peierls phase with 
'\/27r0^ = Tin, where n is an integer. Here, the Kondo 
effect is destroyed due to the spin gap induced by the 
spin-Peierls ordering in the spin chain. As a result, the 
Fermi momentum is kp. 



D. A coexisting Kondo and spin-Peierls phase 

An interesting possibility is described by the RG fixed 
point with yK = ysp = +oo. This is the regime in 
which the system is simultaneously in the Kondo-dimer 
and spin-Peierls phases. Then all fields (/>*''^, 0^''^ at- 
tain nonvanishing expectation values and the spectrum 
of spin excitations becomes gapped. There is still a gap- 
less charge-density mode at 2k*p in this state, though. 

Because of the broken translational symmetry, the co- 
existing phase and the purely spin-Peierls phase have the 
same Fermi volume (length). They are, however, distinct 
phases due to the absence (presence) of the Kondo-dimer 
order parameter and the spin gap in the spin-Peierls (co- 
existence) phase. 



If we extend our study towards Ising anisotropy of the 
spin chain (see Appendix IX)). we will have to add Ising 
phase to the total picture as well. However, this phase is 
of little interest in the context of our problem. 



V. QUANTUM PHASE TRANSITIONS 

We now consider in some detail the transitions between 
the phases. For simplicity, we will again set Kg = 1. 



A. Transition from the decoupled Luttinger liquid 
to the Kondo phase 

First, we will study the transition from the Luttinger- 
liquid phase to the Kondo-dimer one. For ysp < x, 
X — Kr — 1, the backscattering interaction ygp always 
renormalizes to zero. There is, however, a separatrix be- 
tween the regions where is relevant and where it is ir- 
relevant (Fig. III). By substituting Eq. into Eq. I)17b|) . 
we find that the equation of the separatrix is 



sen ^1 ' 

Vi/ = — In , 

2 V2a; 



(1:) + ^^ 



(22) 



where x*"^ — Xq — Vsp- Thus, the phase boundary cor- 
responding to this transition is determined by the initial 
conditions for isotropic Kondo interaction y± ~ Vz — Uk 
of the form — ifx' {xq). An approximate equation 
of the phase boundary is yx — x*"^ /Au for x* <C u, in 
agreement with the result obtained at the end of Sec. lIIII 
When Kondo interaction is stronger than y^^, the 
flow of the constants y±^ and yz eventually becomes rel- 
evant and X approaches zero. Suppose that initially 
Uk — Uk^ (xq) + t, where t <^ x* . Then the flow of x 
is determined by the following differential equation: 



dx 



In 



X 



+4:Ut In 



(23) 



By expanding it about x — x and assuming that x* <^ u, 
we derive 



dx 



2t 



,*3 



„*2 



{x- 



(24) 



The magnitude of the spin gap is determined by the value 
of the parameter I when becomes of order of unity, 
A ~ exp(— /*). Thus, we find that in the Kondo-dimer 
phase near the transition the spin gap is exponentially 
small. 



Ak ^ exp 



(25) 
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Therefore, the transition from Luttinger-liquid phase 
to the Kondo-dimer phase is continuous and, indeed, be- 
longs to the Kosterhtz-Thouless type. However, the ex- 
ponent is not universal due to the dependence on xq. 
The formula above is correct only for t — yx — ifx^ <C a;* 
and it becomes invalid as x* 0, where it is A/f ^ 
exp {—1/ut). 

Now consider the behavior of the correlation functions 
at the transition. We observe that the correlation func- 
tions (s(x,t) • s(0)) and {t{x,t) ■ t(0)) are nonuniversal. 
Indeed, the uniform part of the former decays as r"'''" and 
the uniform part of the latter is r~^^ , where r = x ± vt 
and 



= 2- 



{Vr + Vsf 
1 



Kr 



" 2 — 3 



(26a) 
(26b) 



Here — Vz (a;) on the separatrix, ~ x'^/2u. The 
staggered part of the localized-spin susceptibility decays 
as 7--i/^T^ g^g involves only Or and no (fir- The corre- 
sponding local spin susceptibility, x(w;x = 0), is diver- 
gent. (The divergence at the QCP becomes logarithmic 
in the spin-isotropic limit, which is consistent with an ex- 
trapolation towards the QCP of the result obtained inside 
the Kondo-dimer phase of the spin-isotropic model using 
the form-factor techniquci^i) 

In order to construct a correlation function that is 
characterized by a universal critical exponent at the tran- 
sition, we will introduce the following operator, defined 
in terms of right- and left-moving fields: 



(27) 



The fields i?^ and L% are the spin pieces of the right- 
and left-moving fermions in the IDEG, introduced ear- 
lier, so that R%, L% ~ exp 



This 

operator mixes the spinous of the spin chain and those 
of the conduction electron gas. While the spin does not 
flip in the process of such transformation, the momen- 
tum is not conserved and changes by k*p. Judging by 
the representation of P in the bosonic fields, one could 
formally identify it with the square root of the Kondo 

1 /2 

singlet, P ^ (sj^ ' "Tl) • The susceptibility associated 
with P, xp(2;,t) = (P (a;, r) P (0)), decays as a power 
law r"^"*^/^, therefore, the corresponding local suscepti- 
bility at the transition (where e = 0) has universal be- 
havior 



XP (w; X = 0) - In 



0, 



(28) 



B. Transition between the Kondo, spin-Peierls, 
and coexistence phases 

If we start at the separatrix Eq. (|22|l and begin to 
decrease Kr, we will get deeper and deeper inside the 
Kondo phase. By iterating the RG equations, we found 
that there are two additional regions of the phase dia- 
gram. 

For Kr sufficiently smaller than the separatrix value, 
we find that renormalizes to zero while ysp renor- 
malizes to larger and larger values. This yields the pure 
spin-Peierls phase discussed in Sec. IIVCI 

Between the Kondo phase and the spin-Peierls phase, 
we find a finite region of parameters over which both 
y±^ and ysp renormalize to larger values and reach order 
unity simultaneously. We interpret this as meaning that 
we enter a state in which the Kondo dimers coexist with 
spin-Peierls phase. 

The exact locations of the corresponding phase bound- 
aries are impossible to determine from the RG equations, 
since the latter become invalid when y^^lovysp^l. 
However, we can determine an approximate location of 
the coexistence phase from the condition that both y± 
and ysp reach 0(1) at certain finite value of I. This line 
has the following shape for finite yx'- 



ysp = X 



-^^y^) , x> 0, 



exp 



exp 



^uyK 



Xc < X < 0, 
1/2' 



(29) 



where A is a cutoff. 



where Xc — —2^uyK- The dependence of the spin-gap 
in the spin-Peierls phase near x^ is A^p ~ Vsp'^^^^ ■ 

The same line for finite ysp ends at the point x — ysp, 
at which all three phase boundaries merge and which 
separates spin-Peierls phase from Luttinger liquid when 
yK vanishes. 



C. Phase diagram 

These results specify the phase diagram. It is hard 
to show the diagram in the full parameter space (Jk, 
JsP, Kr), as it would require a three-dimensional graph. 
Instead, we plot two typical cross-sections, corresponding 
to the cut at a finite constant value of Js p (Fig. and 
at a finite constant value of Jk (Fig. . 

One can also easily generalize these phase diagrams 
to the case when Jsp < (see Appendix IXI) . There is a 
complete symmetry with respect to the change of the sign 
of Jsp, except that one will have to replace spin-Peierls 
phase with Ising one on the phase diagrams. 
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FIG. 3: Phase diagram for fixed finite Jsp- The sohd line 
shows the continuous transition and the shaded area shows 
the possible location of the coexistence region. The dashed 
line has been defined in Eq. I|29|l . but the exact shape of the 
area is unknown yet. The lines merge at a single point. The 
brackets label the Fermi momentum of each phase. 




FIG. 4: Phase diagram for fixed finite Jk- The lines and the 
shaded area have the same meaning as in Fig. |3 



VI. CONCLUSIONS 

We have studied a onc-dimcnsional Kondo lattice 
model that consists of a one-dimensional electron gas 
coupled to a spin-1/2 chain through Kondo interaction. 
The phase diagram is quite rich and contains phases with 
either a large Fermi momentum or a small Fermi momen- 
tum. 

When both nearest-neighbor and a significant next- 
nearest-neighbor interactions are present in the isotropic 
spin chain and Kondo interaction is small, the system is 



in the spin-Peierls phase, which is characterized by bro- 
ken translational symmetry due to the presence of spin 
dimers and by small Fermi momentum kp oi conduction 
electrons. 

Increasing the strength of Kondo coupling constant, 
we have found indications that the system enters a co- 
existence phase, containing both the dimers of Kondo 
singlets and the dimers of the spins in the spin chain. 
In this phase the translational symmetry is also broken, 
but the conduction electrons are part of the system with 
large Fermi momentum kp. Additionally, the spin gap 
encompasses conduction electrons as well. 

As we further increase Kondo interaction, the system 
enters a pure Kondo-dimer phase with no translational 
symmetry broken and with large Fermi momentum kp. 
We were unable to determine the precise nature of the 
transitions either between the spin-Peierls phase and the 
coexistence phase or between the coexistence phase and 
the pure Kondo phase, due to the inherent limitations of 
the renormalization-group approach. 

We have also studied the model in which the spin 
chain contains an XY anisotropy. For sufficiently weak 
Kondo interaction, the system is always in a Luttinger- 
liquid state, characterized by small Fermi momentum 
kp and no breaking of translational invariance. Increas- 
ing Kondo interaction triggers a continuous Kosterlitz- 
Thouless transition into the Kondo phase. We have cal- 
culated the dependence of the spin gap on Kondo inter- 
action near the transition and have found that it is ex- 
ponentially small. We have also identified a correlation 
function that logarithmically diverges at zero frequency 
at the transition. 

The possibility of a continuous transition from small 
to large Fermi momenta looks puzzling at first sight. In 
our case, the weight of the Kondo resonance is charac- 
terized by the Kondo-dimer order parameter. The fact 
that the transition is continuous reflects the continuous 
onset of this weight at the transition. However, even in- 
finitesimally small order parameter is still a macroscopic 
quantity. The electron count includes the localized spins 
as soon as the Kondo-dimer order parameter is devel- 
oped, but does not include them as long as we are in the 
Luttinger liquid phase. A jump in the Fermi momentum 
then takes place at the transition point. 

In dimensions higher than one, Fermi-momcntum- 
changing transitions have been discussed in the liter- 
ature. One such transition arises in an extended dy- 
namical mean field treatment of the Kondo latticefiii^S 
between a large-Fermi-momentum paramagnetic metal 
phase and a small-Fermi-momentum antiferromagnetic 
metal phase. Here the transition is continuous and is ac- 
companied by a logarithmically divergent local spin sus- 
ceptibility. The transition is locally quantum critical, in 
the sense that Kondo resonances are destroyed at the 
transition. A related Fermi-momentum-changing transi- 
tion has been discussed in certain large- limit of frus- 
trated Kondo lattice systems, ^'^ between a large-Fermi- 
momentum paramagnetic metal phase and a fractional- 



ized small-Fermi-momentuni paramagnetic phase. The 
asymptotically exact study we have carried out in the 
one-dimensional Kondo lattice model reveals a Fermi- 
momentum-changing transition that bears strong simi- 
larities to the local quantum criticality proposed earlier 
for Kondo lattice systems in dimensions higher than one. 

We close by noting that some of the phases we have dis- 
cussed for the purely one-dimensional Kondo lattice may 
be of direct experimental significance. For instance, the 
organic material (perylene)2[Pt(S2C2(CN)2)2] is believed 
to be a realization of a quasi-one-dimensional Kondo 
lattice^ the system also displays spin-Peierls ordering. 
It would be very interesting to study the Fermi-surface 
properties of this material, as well as the phase transi- 
tions of the system by tuning, say, pressure. 
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APPENDIX A: XY AND ISING ANISOTROPY. 

In our study of the Lagrangian Eq. ()l()(l . we were 
assuming that the spin chain had XY anisotropy. In 
this appendix we will generalize our results for Ising 
anisotropy as well. 

First of all, let us determine the boundary between 
the regions with XY and Ising anisotropy, i.e., the area 
with SU (2) symmetry. In the full phase diagram in the 
space of Jsp, Jk, and this area would be a surface. 
However, if we look at the cross-section for fixed Jki this 
area reduces to a line (Fig.O. 

For Jk = 0, this line (a) actually consists of two pieces. 
One (ai) that separates regions I He and IVb has a stable 
fixed point at the origin, corresponding to the isotropic 
Luttinger liquid. One (02) that separates regions / and 
Ila has a stable fixed point at infinity, corresponding 
to spin-Peierls phase. The regions / and // (both with 
Jgp > 0) flow towards spin-Peierls flxed point, the re- 
gion /// flows towards Ising fixed point on the line 72 
at infinity, which corresponds to Neel order with broken 
Ising symmetry, and the region IV flows towards the line 
of Luttinger-liquid fixed points (32- This is a standard 
Kosterlitz-Thouless picture. 

As Jk ^ 00, the line with SU{2) symmetry rotates 
so that it overlaps with Jgp axis. Now in all regions 
the fiow is towards the line of Kondo-phase fixed points 
along Jsp axis. For finite values of Jk, spin-Peierls and 
Kondo fixed points still possess SU{2) symmetry, while 
Luttinger-liquid and Ising fixed points break it. 

The continuous transition from the Luttinger-liquid to 
the Kondo phase is always in the XY region and it ex- 
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FIG. 5: XY and Ising anisotropy (cross-section at Jk ~ 0). 
The line with SU{2) symmetry is a, the regions II and IV 
have XY anisotropy, and the regions I and III have Ising 
anisotropy. For Jk 00, the line with SU(2) symmetry 
becomes Jsp axis. 



ists for any sign of Jsp, but only for Kt > 1, which 
is why most of the calculations in this paper have been 
performed assuming XY anisotropy. The line segment 
71 that was separating Luttinger liquid from spin-Peierls 
phase for Jk = now splits into two boundaries that 
surround the Kondo phase (Fig.Q)). 

As for the transition from Kondo to spin-Peierls phase, 
its study is more complicated due to the mentioned possi- 
bility of the coexistence region. However, it is clear that 
the finite Jk makes this transition (or the coexistence 
region) extend towards a finite range within region / as 
well. We have determined the location of the coexistence 
region for Ising anisotropy with Jgp > and K^- < 1, 
which is refiectcd in Eq. H29|) . 

What will happen for the remaining region with Ising 
anisotropy {111)1 Let us observe that the RG equations 
Eq. (|13|l are invariant with respect to the change of the 
sign of Jsp- This means that both of our phase diagrams 
Figs. 13141 should remain the same, except that now Jsp 
will need to be replaced with —Jsp and Ising order will 
take place of spin-Peierls one. In the coexistence region 
(now for Kondo and Ising phases) the electrons not only 
will become spin-gapped, but also will break SU (2) in- 
variance so that Ising order will effectively expand over 
the entire system. 
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APPENDIX B: SMALL FERMI MOMENTUM IN 
THE DECOUPLED LUTTINGER LIQUID 

In this appendix, we consider the decoupled Luttinger 
hquid phase in which the transverse component of the 
Kondo couphng [Jj^) is irrelevant. We wish to establish 
that the Fermi momentum is kp- This statement is ob- 
viously true if the longitudinal component of the Kondo 
coupling is absent, in which case the conduction 

electrons would be completely decoupled from the spin 
chains when we reach the fixed points. Below, we show 
that it remains valid even for a finite J^^-coupling. 

The J^j^ term describes the forward-scattering inter- 
action between the z— components of the spins of the 
conduction electrons and those of the spin-chain [cf. 
Eq. Hllc|l ]. It is marginal in the RG sense when = Q 
[cf. Eq. (|13bp ]. In addition, both Kr and Kg are also 
marginal since, in the decoupled Luttinger liquid phase, 
ysp = as well. We can then introduce the following 
diagonalization (mentioned already in Section llll(l : 



0s = as(l)s + bs(j)T, 

(j)'^ = ar(t>s + br(j)T, 
v'g = CsVs + dgVr, 
v'^ — CtVs + drVr- 



(Bla) 
(Bib) 
(Blc) 
(Bid) 



The transformation coefficients can be straightforwardly 
derived for arbitrary values of Jj^ . For simplicity, how- 
ever, we will write down the expressions that are valid 
only up to the second order in J' 



K ■ 



a, = 1 



(i;2_^2)2— V27r 



VsK. ( J. 



Cc 1 



vl V27r 



2(«2_,2)2 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



The expressions for a,-, Cr, and d-r have the same 
forms, except that the subscripts s and t are exchanged. 

The single electron Green's function, Ga{x,t) = 
— (TtCa{x,t)cl^ {0,0)'^, can now be calculated using the 
bosonization form: c^ix) = X]r=± '^'"^ (^)' where r — ± 



labels the right/left moving parts and 
1 



^^rk^xp^^ exp {-ir<i>raix)} . (B3) 



Here, ^raix) = \fT^\4>a + rOa\- For definiteness, we con- 
sider the Green's function of a right-moving electron, 

G+,{x,t) = -(rtc+,(2;,t)4,(0,0)). Using - 

(l/\/2)<i>c,+ -l-(cr/-\/2)$s,+ and the properties of the Klein 
factors described in the main text, we can see that 
G-(-(x,t) factorizes into a charge part and a spin part. 



G+,(a;,i) = --^e^'^--"Ge,+ (x,t)Gs,+,(a:,i), 

ZTTfl 



G,,+(a;,i) = (Ttexp<i--^$,,+ (a;,i) 



xexp|-^<i>,,+ (0,0)!>), (B4) 



Gs,+aix,t) = ( Ttexp \ - — <I>s,+ (a;,i) 



X exp <j -^$s,+ (0,0) 



The charge part, Gc.\{x,t)^ is determined by the IDEG 
alone. Because of the J^^ couphng, the spin part, 
Gs,+a{xi t), does involve the spins of both the IDEG and 
the spin-chain. This coupling, however, can be han- 
dled by going to the diagonalized basis (f)'^ and 4>'^ intro- 
duced in Eq. (jBip . Since the diagonalization affects only 
the q spin component of the spin chain, going to the 
primed basis does not introduce any oscillatory factor in 
the spatial dependence of Gs,+a{x,t). As a result, the 
single-electron Green's function G+o-(a;, t) is a product of 
gikpx g^j^j factors which decay algebraically in x — Vct 
and X — Vst, respectively. The Fermi momentum, hence, 
is kp. 

The key to the above reasoning is that only the forward 
scattering component of the longitudinal Kondo interac- 
tion is marginal. All other components, including those 
that involve either the q = ±7r mode of the spin chain 
or the q = ±2kp spin mode of the IDEG [cf. Eqs. P7|l ] 
are irrelevant. In particular, the irrelevant nature of the 
interactions involving the staggered moment of the spin 
chain leaves no room for the single electron Green's func- 
tion to contain any Fermi momentum other than fcp. 
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